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ABSTRACT 

Using semi-analytic modeling techniques, we calculate the luminosity function of 
galaxy populations residing in cold dark matter halos of different mass. We pay par- 
ticular attention to the influence of the reionization of the Universe on the number of 
faint galaxies and to the effects of dynamical friction and tidal limitation of satellites 
on the number of bright galaxies. We find substantial differences in the shapes of the 
galaxy luminosity functions in halos of different mass which reflect generic features of 
the cold dark matter model of galaxy formation and thus offer the opportunity to test 
it. We then consider how the individual halo luminosity functions combine together to 
produce the global luminosity function. Surprisingly, the global function ends up hav- 
ing a shallower faint end slope than those of the constituent halo luminosity functions. 
We compare our model predictions with the limited datasets compiled by Trentham 
& Hodgkin (2002). We find good agreement with the luminosity functions measured 
in the Virgo and Coma clusters but significant disagreement with the luminosity func- 
tions measured in the Local Group and Ursa Minor cluster. We speculate on possible 
inadequacies in our modeling and in the existing observational samples. The luminos- 
ity functions of galaxies in groups and clusters identified in the 2dF and SDSS galaxy 
redshift surveys offer the prospect of testing galaxy formation models in detail. 



1 INTRODUCTION 

Understanding the galaxy luminosity function has been a 
goal of galaxy formation theory for several decades (e.g. 
White & Rees 1978; Cole 1991; White & Frenk 1991). 
A particularly interesting question is whether the luminos- 
ity function is universal or whether it depends on environ- 
mental factors such as the mass of the dark halo that hosts 
a particular galaxy population. Considerable attention has 
been paid to the faint end of the luminosity function which 
has a much flatter slope than the low mass end of the halo 
mass function predicted in cold dark matter (CDM) models 
of galaxy formation (e.g. Norberg et al. 2002). 

The early work of White & Rees (1978) showed that 
the number of faint galaxies must have been strongly af- 
fected by feedback processes that prevented most of the gas 
from cooling in small halos at early times. Some likely feed- 
back mechanisms, such as the injection of energy into the 
interstellar medium in the course of stellar evolution, de- 
pend on the internal properties of the galaxy and so their 
effects may be expected to be independent of the large-scale 
environment. A number of observational studies, such as a 
recent analysis of the 2dF galaxy redshift survey (De Propris 
et al. 2002), indeed find no significant difference between the 
luminosity functions of galaxies in rich clusters and in the 



field. Other studies, however, have found the opposite. For 
example, Phillipps & Shanks (1987) concluded that galaxies 
in rich clusters have luminosity functions with considerably 
steeper faint ends than galaxies in the field. More recently, 
Trentham & Hodgkin (2002) have claimed that the faint 
end of the galaxy luminosity function varies systematically 
with environment, increasing in slope from small, diffuse sys- 
tems like the Local Group, to massive, dense systems like the 
Coma cluster. 

In the cold dark matter model of galaxy formation, dark 
matter halos retain considerable substructure after they col- 
lapse and virialize (e.g. Klypin et al. 1999; Moore et al. 
1999) and some of these subhalos are associated with sites 
of galaxy formation. The mass function of subhalos appears 
to be relatively independent of the mass of the parent halo. 
Thus, trends such as those inferred by Trentham & Hodgkin 
(2002) would require processes that either preferentially sup- 
press the formation of dwarf galaxies in low mass systems, 
or destroy them after they form. An effective mechanism 
for suppressing the formation of small galaxies is the re- 
heating of the intergalactic medium (IGM) caused by the 
reionization of the Universe at a redshift z ~ 6. Tully et al. 
(2002) have argued that this process could introduce an en- 
vironmental dependence in the galaxy luminosity function 
on the grounds that a higher fraction of the low-mass halos 
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that formed before reionization (when dwarf galaxy forma- 
tion proceeded unimpeded by photoionization suppression) 
ended up in clusters today than in less massive systems. 

The effect of reionization on the formation of galax- 
ies has been the subject of several recent studies (Bullock, 
Kravtsov & Weinberg 2000; Benson et al. 2002a; Somerville 
2002), aimed mostly at investigating the discrepancy be- 
tween the large number of subhalos found in N-body simu- 
lations of galactic CDM halos and the small number of satel- 
lite galaxies observed in the Local Group. In this paper, we 
employ the CDM model of Benson et al. (2002a) to calcu- 
late the luminosity function of galaxy populations residing 
in dark matter halos of different mass. We find that there 
are significant differences in these luminosity functions and 
we then explore how they combine together to build up the 
global luminosity function, with particular emphasis on the 
faint end slope. A partial study of luminosity functions in 
halos of different mass using a semi-analytic model of galaxy 
formation was carried out by Diaferio et al. (1999). 

To calculate galaxy luminosity functions in halos of dif- 
ferent mass correctly, it is important to include tidal effects 
on satellite galaxies, a potential galaxy destruction mecha- 
nism. Our model treats these effects in considerably more 
detail than previous models of galaxy formation. We find 
that tidal effects are important in limiting the formation of 
massive galaxies at the centre of rich clusters. 

In this paper, we compare the results of our calcula- 
tions to the data of Trentham & Hodgkin (2002) and assess 
whether feedback from reionization is a viable explanation 
of the trend claimed by these authors. The existing dataset 
is small, but forthcoming results from the 2dF and Sloan 
galaxy surveys will enable much more extensive comparisons 
with the theory. 

The remainder of this paper is arranged as follows. In 
§2 we briefly outline our model of galaxy formation, in §3 
we present results for the environmental dependence of the 
luminosity function and in §4 we compare our model with 
the available observational data. Finally, in §5 we present 
our conclusions. We present, in an Appendix, several simple 
models of photoionization suppression to elucidate how this 
mechanism works. 



2 MODEL 

We employ the semi-analytic model of galaxy formation 
described in detail by Cole et al. (2000) and Benson et 
al. (2002a; hereafter Paper I) to compute the properties of 
galaxies in a range of environments at z = 0. The reader 
is referred to those papers for a complete and detailed de- 
scription of the model. Briefly, the hierarchical formation of 
dark matter halos is calculated using the extended Prcss- 
Schechter formalism (Press & Schechter 1974; Bower 1991; 
Bond et al. 1991). The formation of galaxies in the resulting 
dark matter halo merger trees is followed by means of simple, 
physically motivated models of gas cooling, star formation 
and galaxy merging. Recent work has demonstrated that at 
least in so far as gas cooling is concerned these simplified 
calculations are in excellent agreement with the results of 
N-body/hydrodynamical simulations (Benson et al. 2001; 



Yoshida et al. 2002; Helly et al. 2002). Applying a stel- 
lar population synthesis model gives galaxy luminosities in 
different passbands. The model includes a prescription for 
supernovae feedback which drives gas out of galaxies at a 
rate proportional to the current star formation rate, with 
a constant of proportionality that is larger for less strongly 
bound systems. This negative feedback flattens the faint end 
slope of the luminosity function but, on its own, its effect 
is not strong enough to account for the measured function 
(Paper I). Benson et al. (2002a) developed, and incorpo- 
rated into this model, a detailed, self-consistent treatment 
of photoionization suppression of galaxy formation (here- 
after abbreviated as PhS). By causing further suppression 
of faint galaxy formation, this mechanism brought the model 
luminosity function into excellent agreement with observa- 
tions. In Paper I, we also included a much more detailed 
treatment than was previously possible of the evolution of 
satellite galaxies as they orbit within larger halos experienc- 
ing tidal forces and gravitational shock heating. (Hereafter, 
we abbreviate "tidal limitation as TiL.) The associated mass 
loss also affects the shape of the luminosity function, and so 
we will explore this process in this paper. 

We adopt essentially the same model parameters as in 
Paper I, with the following differences. While the models of 
Paper I used the Press-Schechter halo mass function, we 
choose to adopt the function proposed by Jenkins et al. 
(2001) which gives a better match to the results of N-body 
simulations. Consequently, to produce a reasonable galaxy 
luminosity function, we find it necessary to adjust the pa- 
rameters O b from 0.020 to 0.024, e* (which determines the 
star formation rate) from 0.0050 to 0.0067, and T (which 
affects mass-to-light ratios) from 1.38 to 1.03 (since T af- 
fects the fraction of mass recycled in star formation, we ad- 
just that fraction accordingly). Even with these changes, the 
Jenkins et al. (2001) mass function results in a somewhat 
worse fit to the bright end of the luminosity function than 
was obtained in Paper I. We defer further study of this as- 
pect of the model to future work. We adopt an escaping 
fraction of ionizing photons from galaxies, / csc , of 12%. This 
produces a neutral hydrogen fraction for gas at zero overden- 
sity in the IGM at z = 6 of 5 x 10~ 4 , in agreement with the 
lower limit of 3.4 x 10~ 4 derived by Lidz et al. (2002), while 
causing reionization to occur at as high a redshift as pos- 
sible. This maximizes the associated suppression of galaxy 
formation*. The resulting model is essentially the same as 
the / csc = 10% model presented in Paper I. 



3 THE GALAXY LUMINOSITY FUNCTION 

3.1 Luminosity functions in halos of different mass 

We begin by studying the luminosity functions of galax- 
ies in individual dark matter halos of different mass, and 
then proceed to consider how these combine together to 

* We could increase / eS c slightly further without violating the 
limit of Lidz et al. (2002) but, since the neutral hydrogen fraction 
drops very rapidly in our model at these redshifts (Paper I), it 
makes almost no difference to the suppression of galaxy formation. 
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form the global galaxy luminosity function. Using the model 
described in §2, we simulate galaxy formation in a large 
number of halos of masses 10 11 , 10 12 , 10 13 , 10 14 and 
10 15 h~ 1 MQ^ . This spans halo masses appropriate to sys- 
tems ranging from the Local Group to rich galaxy clus- 
ters. (The 10 11 ft -1 Mq bin corresponds to halos less mas- 
sive than that of the Local Group but we include it since 
it shows the most dramatic effects.) Merger trees are con- 
structed with a resolution of 10~° of the final halo mass or 
10 8 /i _1 Mq, whichever is smaller. This is sufficient to follow 
accurately the formation and evolution of galaxies brighter 
than Mb — 5 log ft = — 10 except in the most massive clus- 
ters (10 15 /i^ 1 Mq) for which resolution begins to affect our 
results faintwards of Mb — 5 log ft = —12; for these clus- 
ters, we show results only down to Mb — 5 log h = —13. 
We perform calculations for our "standard model" (i.e. in- 
cluding the effects of PhS and TiL of satellite galaxies), and 
repeat them switching off PhS and/or TiL (resulting in a 
total of four different sets of results) , in order to assess the 
effects of each physical mechanism separately. We simulate 
a sufficiently large number of halos of each mass so that the 
average luminosity functions are reasonably smooth. 

The solid lines in Figure 1 show the luminosity func- 
tions of galaxies in our simulated halos. The basic form of 
the standard model luminosity function is similar in all mass 
halos. It has two components: a near power-law distribution 
at faint magnitudes which is made up of satellite galaxies 
(i.e. those galaxies that do not reside at the centre of the fi- 
nal halo), and a strongly peaked distribution at bright mag- 
nitudes beyond which no more galaxies are found and which 
is made up of the central galaxy in each final halo. As we 
consider halos of increasing mass, both components shift to 
ever brighter luminosities, but the peaked component de- 
creases in amplitude relative to the power-law component, 
indicating that, in low mass halos, a single galaxy tends to 
dominate, while in clusters the light is more equally shared 
among many galaxies. 

Observed luminosity functions, 4>{L) — dn/dL (where 
n is the number of galaxies per unit volume and L is lumi- 
nosity), are typically described by the Schechter function, 

-*.(£)"-(-£)£ 111 

where (/>,, L* and a are parameters. At magnitudes faint 
compared to the characteristic luminosity, L», has the sim- 
ple form 4>(L) oc L a . Following normal practice, we will re- 
fer to a as "the faint end slope." Although the luminosity 
functions in Fig. 1 are not well fit by the Schechter form, 
we can still obtain a formal value of a by fitting to the 
data in the first few bins plotted. We find slopes in the 
range —1.57 < a < —1.48, except in the lowest mass ha- 
los (10 11 h~ 1 M & ) for which a much flatter slope, a = —1.06, 
is found (although it is possible that a steeper slope would 
be obtained also in this case if we probed to even fainter 
luminosities). There is little evidence for a strong variation 
of slope with halo mass, except perhaps for a weak (and 

t Throughout this paper we write Hubblc's constant as Hq = 
100/ikm s" 1 Mpc -1 . 



rather insignificant) trend for flatter slopes in higher mass 
clusters. However, it is important to note that the slope is 
not constant with galaxy magnitude. This is clearly seen, for 
example, in the 10 12 ft -1 Mq halos for which, in the region 
— 18 < Mb — 5 log ft < — 15, the slope is close to a = — 1. 

Switching off TiL in our model, and reverting to the sim- 
pler calculation of dynamical friction of Cole et al. (2000), 
results in the dotted curves of Fig. 1. Clearly, TiL has lit- 
tle influence on the faint end of the luminosity function, 
although it does slightly reduce the luminosities of already 
faint galaxies (as is particularly noticeable for lO 12 ft _1 M0 
halos). A more interesting effect of TiL is apparent in the 
massive clusters (10 14 and 10 15 ft _1 M© halos). Here, switch- 
ing off TiL results in a new population of very luminous 
galaxies, visible as a "bump" in the luminosity function 
at bright magnitudes. These highly luminous objects form 
through multiple mergers at the centres of cluster halos. 
When TiL is included, the merger rate of massive galax- 
ies in clusters is significantly reduced, and none of these 
highly luminous objects form. As shown in the left-hand 
panel of Fig. 2, cluster central galaxies are about an order 
of magnitude less massive when our accurate calculation of 
dynamical friction is used. Note that the inclusion of TiL 
makes much less difference in 10 12 ft _1 M Q halos, due to the 
very different shape of the luminosity function in these sys- 
tems (i.e. there is relatively little mass in satellites to merge 
into the central object). This difference between the merger 
rates in our model and that of Cole et al. (2000) cannot 
be fully reproduced by simply lengthening the dynamical 
friction timescale in the Cole et al. (2000) prescription by 
some factor (even if this factor is a function of halo mass) . 
The shape of the distribution of merger times is significantly 
different between the models and, in our calculation, varies 
with halo mass. 

An effect similar to that shown by the comparison of 
the solid and dotted lines in Fig. 1 was noticed by Springel 
ct al. (2001) in cluster luminosity functions calculated using 
a different semi-analytic model. Springel et al. (2001) com- 
puted luminosity functions adopting both simple estimates 
of galaxy merger rates (similar to those assumed in com- 
puting the dotted curves in Fig. 1) and merger rates taken 
directly from a high resolution simulation of dark matter. 
They found that the simple scheme produced highly lumi- 
nous galaxies through merging that were not present in the 
N-body calculation. Based on the number of massive merg- 
ers quoted by Springel et al. (2001), we can conclude that a 
similar reduction in central galaxy mass occurred when our 
realistic merger rate calculation was employed. As suggested 
by Springel et al. (2001), the cause of this difference is that 
the simple scheme underestimates the merging timescale of 
rather massive subhalos in clusters, leading to many massive 
galaxies merging with the central object in the final halo. In 
our more detailed calculation, as in the N-body simulation, 
these large subhalos loose a significant amount of their mass 
as they orbit within the cluster, leading to a reduced dy- 
namical friction force and hence to a reduced merger rate. 

The right-hand panel of Fig. 2 shows the mass func- 
tions of surviving substructure halos in clusters at z = 0. 
The dotted histogram shows the model with no TiL, while 
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Figure 1. B-band galaxy luminosity functions in halos of different mass. Each panel shows the mean predicted model luminosity function 
in an ensemble of dark matter halos of mass as given in each panel (ranging from small galaxies to clusters). Specifically, we show the 
mean number of galaxies per magnitude, per halo. Solid lines include photoionization suppression (PhS) and tidal limitation (TiL), 
dotted lines have PhS only, short-dashed lines have TiL only, and long-dashed lines have neither. Open circles show observed luminosity 
functions from the compilation of Trcntham & Hodgkin (2002), as indicated in the legend of each panel. These have been normalized 
arbitrarily to permit easier comparison of their shapes with the models. Filled squares in the 10 12 h~ 1 Mq panel show the luminosity 
function compiled by Benson et al. (2002b), which includes only galaxies classed as lying within the virial radii of the Milky Way's or 
M3f 's dark halos. The comparison between the model and the observational data is discussed in §4 
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the dashed histogram corresponds to the model with TiL 
included, but the mass of each substructure is plotted as it 
was before it experienced any tidal mass loss. Clearly, in- 
clusion of TiL has resulted in fewer mergers since there are 
more surviving substructures of a given initial mass. The 
solid histogram gives the distribution of substructure masses 
after tidal mass loss. Comparing the dashed and solid his- 
tograms indicates that halos typically lose between 50% to 
70% of their mass (with higher mass halos losing less than 
lower mass halos). The dynamical friction time scales as the 
mass of the satellite halo and is therefore increased by fac- 
tors of 2 to 3 when TiL is included, greatly reducing the 
galaxy merger rates. 

Switching off PhS (short-dashed lines in Fig. 1) has 
a more dramatic effect on the luminosity functions than 
switching off TiL. Firstly, the faint end slope steepens in 
almost all systems to values in the small range —1.71 < a < 
— 1.64. (The exception are the lowest mass halos which have 
a — —1.38 because of the particularly strong effects of super- 
novae feedback in these halos.) The increase in the number 
of faint galaxies is more pronounced in lower mass systems. 
In high-mass systems there is a weak trend for the bright 
end of the luminosity function to be shifted slightly faint- 
wards when PhS is ignored. In this case, more gas is locked 
up in the fainter galaxies, leaving less gas to form massive, 
bright galaxies. In the lowest mass systems the opposite ef- 
fect occurs and the central galaxies become brighter when 
PhS is ignored. In these 10 11 ft -1 Mq halos, PhS suppresses 
the amount of gas that is able to accrete into the halo, thus 
reducing the luminosity of the associated galaxy. The way 
in which various processes associated with photoionisation 
shape the luminosity function in halos of different mass, par- 
ticularly its faint end, is discussed in detail in the Appendix. 
There we also discuss the effect of various approximations 
in the treatment of photoionisation. 

Long-dashed lines show results when both TiL and PhS 
are left out of our calculations. These are essentially accumu- 
lations of the effects seen when each process was switched off 
individually, and result in a further steepening of the faint 
end slopes. In these models, only feedback from supernovae 
affects the faint end of the luminosity function in a manner 
which is essentially independent of halo mass. 

3.2 The global luminosity function 

We now explore how the individual halo luminosity func- 
tions presented in Section 3.1 combine to produce the global 
luminosity function. An initially surprising aspect of the in- 
dividual halo luminosity functions is that their faint end 
slopes are, in all cases, steeper than that of the global lumi- 
nosity function. This, as we showed in Paper I, has a slope 
of around —1.2 both in our models and in the real universe. 
How does PhS contrive to produce a flatter faint end slope 
for the total population than for any of the constituent in- 
dividual halos? To understand this behaviour, we must con- 
sider the contribution of different mass halos to the global 
luminosity function. 

The upper panels of Fig. 3 show the contributions to the 
global luminosity function from satellite and central galax- 



ies (left and right-hand panels respectively) residing in ha- 
los of different mass. The slope of the combined satellite 
luminosity function reflects the power-law slope of the in- 
dividual halo luminosity functions and is dominated over a 
wide range of luminosities by halos of 10 13 -10 14 /i _1 M . The 
luminosity function of central galaxies, on the other hand, 
has a much flatter faint end slope, which is determined by 
the way in which central galaxy luminosity scales with halo 
mass. In our models, this scaling depends primarily on the 
combined effects of supernovae and PhS feedback. Since the 
central galaxy luminosity functions of individual halos are 
strongly peaked around a particular luminosity, the global 
central galaxy luminosity function is dominated by halos in 
a narrow range of mass at each luminosity. (The peak for 
galaxies in 10 10 ' 5 to 10 11 ' 5 h^ 1 Mq is much broader, since 
these mass scales are affected by PhS, causing strong sup- 
pression of galaxy formation in the lower mass halos in the 
range. Figure 1 shows much more clearly the sharply peaked 
luminosity distribution of central galaxies in halos of fixed 
mass.) 

The lower left-hand panel of Fig. 3 shows how the satel- 
lite and central galaxy luminosity functions combine to pro- 
duce the global luminosity function (which is compared to 
the recent observational determination from the 2dF galaxy 
redshift survey of Norberg et al. 2002). Brightwards of 
Mb — 51og/i £s —17, central galaxies dominate, while faint- 
wards of this, the satellite contribution gradually takes over. 
The faint end slope of the global luminosity function ends up 
being intermediate between that of the satellite and central 
galaxy luminosity functions. 

Finally, in the lower right-hand panel of Fig. 3 we show 
the contribution to the global luminosity function from halos 
of different mass. Clearly, the faint end slope of the global 
luminosity function is determined in part by the relation 
between central galaxy luminosity and halo mass, not just 
by the faint end slopes of the individual halo luminosity 
functions. This explains why the global galaxy luminosity 
function ends up having a flatter faint end than the individ- 
ual halo luminosity functions. It is interesting to note that 
both the total and cluster (i.e. 10 1E 7i _1 Mq halo) luminos- 
ity functions are reasonably well fit by Schechter functions, 
albeit with slightly different parameters (a = —1.31/ — 1.45 
and M* — 5 log /i = -20.12/ - 19.60 for the total/cluster 
luminosity function; note that the correlation between the 
parameters a and M* should be borne in mind when com- 
paring these values). 



4 COMPARISON WITH OBSERVATIONS 

Our model predicts marked differences in the luminosity 
functions of galaxy populations residing in halos of differ- 
ent mass. Although the detailed form of these predictions 
depends on the details of our galaxy formation model, the 
gross differences seen in Fig. 1 are generic predictions that 
bear directly on basic features of the model such as hier- 
archical clustering from CDM initial conditions, gas cool- 
ing, photoionization, feedback, etc. Thus, in principle, in- 
dividual halo luminosity functions provide a strong test of 
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Figure 2. Left-hand panel: The distribution of baryonic masses of galaxies residing at the centres of 10 15 /i _1 Mq clusters at z = 0. The 
solid histogram shows results from our standard model which includes TiL while the dotted histogram corresponds to a model with no 
TiL. Right-hand panel: The mass function of surviving substructure halos in 10 15 /i _1 Mq clusters at z = 0. The dotted histogram is 
from a model with no TiL, while dashed and solid histograms are from the model with TiL and give the masses of halos before and after 
tidal mass loss respectively. 



the CDM galaxy formation paradigm. Unfortunately, imple- 
menting such a test in practice is difficult because of the ob- 
servational complications inherent in identifying galaxies at- 
tached to a particular dark matter halo. Gravitational lens- 
ing is a powerful and promising tool for detecting dark mat- 
ter halos directly, but this technique is still in early stages 
of development (e.g. Mellier 2002). A less direct, but still 
useful approach, consists of finding groups and clusters in 
large redshift catalogues such as the 2dF and SDSS surveys. 
When interpreted with the aid of cosmological simulations, 
it is possible to go some way towards identifying galaxy pop- 
ulations likely to be associated with single halos (and their 
subhalos) of a given mass (Eke et al. 2002) . This approach 
will yield interesting data for comparison with our model 
predictions in the near future. 

In the interim, Trentham & Hodgkin (2002) have pro- 
vided a limited compilation of observational data. They give 
estimates of the B-band luminosity functions of galaxies in 
the Local Group and in the Ursa Major, Virgo and Coma 
clusters. In our model, these systems are expected to re- 
side in halos of mass similar to those we have simulated, 
i.e. 10 12 , 10 13 , 10 14 and lO 15 ft _1 M respectively (Benson et 
al. 2002a; Trentham & Hodgkin 2002; Schindler, Binggeli 
& Bohringer 1999; Geller, Diaferio & Kurtz 1999). We 
show Trentham & Hodgkin's (2002) compilation in Fig. 1 
as open circles. The vertical normalization of the datasets 
have been adjusted arbitrarily to permit an easier compar- 
ison of the shapes and faint end slopes of the model and 
observed luminosity functions. For the two most massive 
systems in the figure, our standard model is in reasonably 
good agreement with the observed luminosity function, in- 



cluding the faint end slope, but for Ursa Major and the Lo- 
cal Group, the observations indicate a much shallower slope 
than is produced by our model. The same discrepancy was 
noted by Benson et al. (2002b) for the Local Group, and 
is also apparent in the predicted Local Group luminosity 
function of Somerville (2002) (which has a ~ —1.5 over the 
range —20 < My — 5 log hro < —10, where /170 is the Hubble 
constant in units of 70km/s/Mpc, for the case of reioniza- 
tion at z = 8 which is the closest to our own calculations), 
suggesting that this is a rather robust prediction of PhS. 

In making the comparison between theory and observa- 
tions in Fig. 1, one should keep in mind the possibility that 
the two may not correspond to quite the same quantity. 
The theoretical predictions pertain to luminosity functions 
of galaxy populations residing in individual dark matter ha- 
los of a given mass. However, there is no guarantee that 
the samples selected by Trentham & Hodgkin (2002) do in- 
deed come from individual dark matter halos. For example, 
these samples might be contaminated by contributions from 
secondary halos which happen to be nearby the dominant 
halo. This could be particularly important for the Ursa Ma- 
jor cluster which is a rather diffuse object and is too large, 
given its mass, to be associated with a single, virialized halo 
(according to the standard theoretical definitions). There 
are ambiguities regarding the Local Group data as well. As 
argued by Benson et al. (2002b), many dwarfs nominally as- 
sociated with the Local Group lie beyond the region which, 
according to our model, contains the (distinct, but similar 
mass) halos of the Milky Way and M31. In Fig. 1, we show 
as filled squares the luminosity function of the galaxies that 
Benson et al. (2002b) would class as lying within the virial 
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Figure 3. B-band galaxy luminosity functions. Upper left-hand panel: The luminosity functions of satellite galaxies (i.e. those which 
are not the central galaxy of their host halo). The heavy solid line shows the luminosity functions summed over all halo masses, while 
thin lines show the contributions from halos in different mass ranges as indicated by the key in the lower right panel. Upper right-hand 
panel: The luminosity functions of central galaxies. Line types are as in the upper left-hand panel. Lower left-hand panel: A comparison 
of the satellite (dotted line) and central galaxy (thin solid line) luminosity functions. Also show is the total luminosity function (heavy 
solid line) and the observational determination of the total luminosity function from Norberg et al. (2002). The contributions from all 
types of galaxy (i.e. satellites and centrals) in halos of different mass to the total luminosirty function. Line types are as in the upper 
left-hand panel. 



radius of the Milky Way or M31. This luminosity function 
also shows a fairly flat faint end. Nevertheless, selection ef- 
fects such as these must be accounted for in a more thorough 
comparison of theory and observations. 

The results of Fig. 1 are shown in a different way 
in Fig. 4 which more clearly shows the variation in the 
shapes of the luminosity functions at different luminosi- 
ties. Here, we plot the effective luminosity function slope, 
a e g = L(d 2 N/dL 2 )/(dN/dL), as a function of absolute 



magnitude for all the systems we have simulated. (For clar- 
ity, we show results only for our standard model and for a 
model with both PhS and TiL switched off.) The luminosity 
functions with no PhS or TiL in 10 13 and 1O 14 /i _1 M ha- 
los show an extended region where a e ff is almost constant, 
indicating a power-law luminosity function. With PhS and 
TiL switched on, the slope is shallower, but the luminosity 
function is no longer a power-law over an appreciable range 
of magnitudes. For the other mass systems, the luminosity 
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function is nowhere well described by a power-law and, in 
every case, the inclusion of PhS makes the slope shallower, 
the effect being larger in the lower mass systems. 

Trentham & Hodgkin (2002) characterized variations in 
luminosity function shapes using the "dwarf-to-giant ratio" 
(i.e. the number of galaxies in some range of faint magni- 
tudes relative to the number at brighter magnitudes). We 
compute a similar ratio for our model with PhS and TiL 
and present the results in Table 1. (We use a slightly differ- 
ent magnitude range for the dwarfs because, for our most 
massive halos, our model is not complete over the range of 
magnitudes used by Trentham & Hodgkin.) 

For the two most massive systems, our model results are 
in reasonably good accord with the observational determi- 
nations, but for the lower mass systems there are significant 
discrepancies. This analysis verifies and quantifies the com- 
parison between the shapes of the predicted and observed 
luminosity functions made in Fig. 1. 



5 DISCUSSION 

Using a semi-analytic model of galaxy formation, we have 
calculated the luminosity function of galaxy populations 
contained in halos of different mass. We find large differ- 
ences in the shapes of these luminosity functions. In smaller 
mass halos (M <;10 13 h~ x M©), the luminosity function has 
a "hump" at bright magnitudes and a roughly power-law 
shape at fainter magnitudes. In more massive systems, the 
shape approaches the familiar Schechter form. These dif- 
ferences reflect the relative contributions of "central" and 
"satellite" galaxies. In the smaller mass halos, the bright 
end of the luminosity function is dominated by a single or 
at most a few bright galaxies and the rest of the popula- 
tion consists predominantly of much fainter satellites. In 
the more massive halos, the light is more evenly distributed 
amongst a large number of galaxies. These differences reflect 
the complex interplay between the processes that establish 
the number of halos and subhalos of different mass (i.e. hi- 
erarchical clustering from CDM initial conditions) and the 
processes that light them up (i.e. gas cooling, stellar evolu- 
tion, feedback, etc). 

The global luminosity function is the weighted sum of 
the luminosity functions of individual halos. Interestingly, 
the faint end slopes of the individual luminosity functions 
are always steeper than the faint end slope of the global lu- 
minosity function. This can be understood in terms of the 
relative contributions to the global luminosity function of 
satellite and central galaxies from individual halos, whose 
luminosity functions have steep and shallow slopes respec- 
tively. The slope of the global function is intermediate be- 
tween these two. 

Our model of galaxy formation represents an advance 
over previous models because of its detailed treatment of two 
processes that are important in establishing the luminosity 
function: the suppression of galaxy formation in small ha- 
los at high redshift (PhS), as a result of the photoionization 
of the IGM by early generations of galaxies and quasars, 
and tidal mass loss from satellite galaxies (TiL). PhS is 



particularly relevant at the faint end and has the net ef- 
fect of reducing the number of faint galaxies, giving rise to 
a flatter luminosity function than would otherwise be the 
case. This outcome is far from obvious. For example, a sim- 
plistic model of PhS in which galaxy formation is assumed 
not to occur in low mass halos below the reionization red- 
shift results, in fact, in a steeper luminosity function (see 
Appendix A). In reality, the effects of PhS depend on the 
mass of the halo relative to the time-dependent "filtering 
mass" (defined as the halo mass which accretes only half 
the mass that it would have accreted in the absence of pho- 
toionization; Gnedin 2000) . Galaxies can form even if their 
halos have mass smaller than the filtering mass, but with 
reduced luminosity. The lower the mass of the halo relative 
to the filtering mass, the fainter is the galaxy that forms. 
The filtering mass is an increasing function of time. Thus, 
although lower mass halos form at higher redshift, the red- 
shift where the filtering mass becomes comparable to their 
mass is also higher, increasing the fraction of them which 
experience strong suppression. These two effects combined 
do tend to produce a flatter luminosity function but the de- 
gree of flattening depends on their detailed balance^ . As we 
have shown previously (Benson et al. 2002a; see also Ben- 
son et al. 2002b), PhS on its own cannot account for the 
shallow faint end slope of the global luminosity function. 
Other processes, such as feedback from energy released by 
stellar winds and supernovae, need to be included in order 
to reduce the number of faint galaxies to the observed level. 

The other process that we have treated in detail, TiL, 
is relevant at the bright end of the luminosity function, par- 
ticularly in massive halos harbouring rich clusters. Mass loss 
from satellite galaxies, as they spiral into the centre of their 
halo, increases the dynamical friction time, so reducing the 
merger rate. This prevents the formation of some of the 
highly luminous and unobserved cluster galaxies that tend 
to appear in models of galaxy formation that do not take 
this process into account. 

Keeping in mind the provisos of §4, we find that our 
model fails to reproduce the steepening of the faint end 
slope of the luminosity function with system richness found 
by (Trentham & Hodgkin 2002). The model agrees quite 
well with the data for the Virgo and Coma clusters, but it 
predicts a steeper faint end than is inferred for the Local 
Group and Ursa Minor. Photoionization is an effective (and 
unavoidable) mechanism for suppressing galaxy formation in 
low mass halos, but it does not produce the kind of "envi- 
ronmental" variation seen by Trentham & Hodgkin (2002). 
At the low mass end, this failure of the model is related 
to the discrepancy that we (Benson et al. 2002b, see also 
Somerville 2002), found previously in our CDM models of 
the Local Group. Photoionization has a dramatic effect on 
the overall number of satellites that survive in galactic ha- 
los and can reduce it to the levels seen in the Local Group. 
However, the satellite luminosity function in the models is 



t A model in which a cut-off in circular velocity rather than mass 
is assumed gives results in quite good agreement with our detailed 
calculations (Benson et al. 2002a). 
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Figure 4. Effective slope of model luminosity functions vs. absolute magnitude for galaxy populations residing in halos of different mass, 
as indicated in the panels. Solid lines correspond to the model that includes both PhS and TiL, while the dashed lines correspond to the 
model that includes neither. 
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Table 1. Dwarf-to-giant ratios for galaxy populations in different environments. The magnitude ranges for dwarfs and giants are specified 
in terms of Mb — 51og/i (to maintain consistency with the rest of this work), and correspond to —16 < Mb < —14 and Mb < — 16 
respectively for h = 0.7. 



Halo mass (r'Mg) JV(-15.23 < Mb — 5 log h < -13.23)/AT(M B - 5 log ft < -15.23) Observed 



10 11 0.24 

10 12 0.40 1.25 ±0.35 

10 13 1.85 0.74 ±0.07 

10 14 1.70 1.47 ±0.18 

10 15 1.68 1.40 ±0.11 



considerably steeper than is observed. This discrepancy con- 
trasts with the success of CDM models in reproducing the 
structural properties of the satellites (Stoehr et al. 2002). 

Taking at the face value the discrepancy between our 
models and the data of Trentham & Hodgkin (2002), it is in- 
teresting to speculate about the possible root causes of the 
problem. Leaving aside alternative theories for the nature 
of the dark matter or drastic changes to the initial power 
spectrum of density perturbations specifically designed to 
reduce the amount of small-scale power (e.g. Colm, Avila- 
Reese & Valenzuela 2000; Kamionkowski & Liddle 2000), 
we focus on our modeling of baryonic processes within the 
context of the CDM cosmogony. This is undoubtedly sim- 
plified. For example, the effects of PhS depend upon both 
the time variation of the filtering mass and on the assumed 
functional form of the suppression, both of which are uncer- 
tain. As discussed in the Appendix, reasonable variations 
in our assumed mass dependence of PhS do not seem ca- 
pable of producing a sufficiently flat faint end slope in our 
models of the Local Group. However, we cannot exclude the 
possibility that the simplifications made in our model are 
inappropriate. For example, we neglect the patchy nature of 
reionization, but if low density regions became ionized first, 
galaxy formation would be suppressed more effectively in 
these regions, altering the temporal evolution of the filter- 
ing mass. Localized PhS triggered, for example, by a nearby 
galaxy or quasar could also change our model predictions. 
These processes are poorly understood theoretically, but ob- 
servational evidence (e.g. Adelberger et al. 2002) suggests 
that they could be important in determining the thermody- 
namic state of the IGM at high redshift. 

At a more basic level, there are approximations in 
our model which might be inappropriate in the regime 
of low mass halos. Thus, although comparisons with N- 
body/hydrodynamical simulations indicate that our esti- 
mates of gas cooling rates are accurate for relatively massive 
galaxies (Benson et al. 2001; Yoshida et al. 2002; Helly et 
al. 2002), it has not been shown that the same is also true 
for the low mass galaxies of interest here. Cooling times in 
these objects can become very short and this might con- 
ceivably render incorrect the assumption of cooling from an 
initial, spherically symmetric quasi-hydrostatic equilibrium. 
Further study of gas cooling in low mass systems is neces- 
sary to address this issue. Similarly, it is far from clear that 
the simple rules for star formation in our model, which are 
empirical parameterizations tuned to match observational 
constraints for relatively bright galaxies at low redshifts, re- 



main valid when extrapolated to low mass galaxies and/or 
high redshifts. Other physical constraints on star formation 
not currently included in our models may also be impor- 
tant for these low mass galaxies (e.g. Verde, Oh & Jimenez 
2002). Finally, our modeling of feedback is also simplified 
and, again, its extrapolation to low masses may be unreal- 
istic. Possible sources of feedback such as strong outflows 
or heating by AGN which are currently neglected may be 
important on these scales. 

In conclusion, we have found substantial differences in 
the luminosity functions of galaxies residing in dark matter 
halos of different mass. These result from the interplay of a 
variety of processes that affect the formation of galaxies of 
different luminosity. Although the exact predictions of our 
model no doubt depend on its details, the gross differences 
between the predicted luminosity functions are likely to be 
generic and to be present in a broad class of cold dark mat- 
ter models of galaxy formation. Comparison with existing 
data gives mixed results: fair agreement with the luminosity 
functions in Virgo and Coma but substantial disagreement 
with the luminosity functions in the Local Group and Ursa 
Minor cluster. The much larger and better controlled sam- 
ples that will be forthcoming from cluster analyses of the 
2dF and SDSS galaxy redshift surveys will allow a more 
comprehensive test of the general class of cold dark matter 
models of galaxy formation. 
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APPENDIX A: PHOTOIONIZATION AND THE 
FAINT END OF THE LUMINOSITY FUNCTION 

The faint end of the luminosity function directly reflects the 
outcome of feedback effects. The two sources of feedback in 
our model, supernova-driven winds and photoionization, are 
both important in determining the number of faint galax- 
ies. As we discussed in Paper I, photoionization inhibits the 
formation of faint galaxies for two reasons: it raises the tem- 
perature and pressure of the IGM, curtailing the ability of 
gas to accrete into halos, and it reduces the rate at which 
gas in halos can cool, inhibiting star formation. The first 
of these processes is the dominant effect, ft can be conve- 
niently described in terms of a "filtering mass," Mf, defined 
as the mass of a halo which accretes only half the baryonic 
mass that it would have accreted in the absence of pho- 
toionization. The mass of gas accreted by a halo of mass M h 
is approximated by the following formula (Gnedin 2000), 
based on the results of gas-dynamical simulations: 

f b M b 



where / D is the universal baryon fraction and Mf is the 
filtering mass. 

There are three different aspects to consider when eval- 
uating the effect of the filtering mass on the faint end slope 
of the luminosity function: 

(i) Galaxy formation is significantly suppressed in halos 
less massive than the filtering mass. 

(ii) The filtering mass is an increasing function of time 
(at least over the range of redshifts of interest here). 

(iii) For halos formed at a given redshift, the suppression 
of galaxy formation is greater the smaller the halo mass is 
relative to the filtering mass. 

The impact of these three factors will depend on the 
distribution of formation redshifts of dark matter halos of 
different mass in different environments^. According to the 
extended Press-Schechter formalism, it is always true (for a 
CDM power spectrum) that the typical formation redshift of 
a halo increases as the mass of the halo decreases, and also as 
the mass of the final host in which the halo resides increases. 
Further understanding of halo formation redshift distribu- 
tions may be gained by considering eqn. 2.15 of Lacey & 
Cole (1993) that gives the mass function of progenitor halos 
at any redshift and which we reproduce in slightly different 
form below: 

(6^ 



1 M hoBt 



■SP) 



dM (Z> " ^ M [a 2 (M) — cr 2 (Mhost)] 3//2 



exp 



(5 



c(0)x2 



2[a 2 (M)-a2(M host )] 



M gas = 



(2V3 _ l)M F /M h ]3 ■ 



(Al) 



where 5^ is the critical value of the linear theory fractional 
overdensity for collapse at redshift z, M is the mass of the 
progenitor halo at that redshift, cr 2 (M) is the variance of 
the linear fractional overdensity in a sphere of mass M, and 
-Mhost is the mass of the halo into which that progenitor has 
been incorporated by the present day. 

Consider galaxies living in relatively low mass halos at 
the present day (in which most of the galaxies that make 
up the faint end of the luminosity function reside). The rel- 
ative numbers of halos of masses Mi and M2(< Mi) at 
some particular redshift depends upon <r 2 (Mi) — ff 2 (Mh ost ) 
and ct 2 (M 2 ) — a 2 (M hos t). As M hos t increases, a 2 (M bos t) de- 
creases for a CDM power spectrum. Thus, for sufficiently 
large Mi, ost , the relative numbers of halos of masses Mi and 
M2 tends to a fixed value. For Mhost comparable to Mi, 
the relative numbers of these halos become a very strong 
function of M hos t, and the number of mass Mi halos is ex- 
ponentially suppressed, as indicated by eqn. (A2) . 

These theoretical expectations are clearly manifest in 
Fig. Al, where we plot the positions of galaxies in the 
halo-mass vs. formation-redshift plane. All of the galaxies 
we consider at 2 = are satellites in a larger host halo 
(of mass Mhost) due to the ranges of luminosity and Mh os t 
we have chosen to consider. Therefore, we plot the forma- 
tion redshift and mass for the halo in which the satellite 



8 By "environment" we mean the mass of the halos into which 
that earlier halo has been incorporated by the present day. 
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formed. We adopt the definition of formation redshift of 
Cole et al. (2000), namely a halo is assumed to be newly 
formed if its mass exceeds twice the mass of each of its 
progenitors at their own formation time, and plot the halo 
mass at the formation redshift (it may have increased after- 
wards due to continued accretion and merging). (The for- 
mation redshifts appear quantized in this figure because of 
the finite timesteps in our calculations which, however, arc 
sufficiently fine so as not to affect the results.) The halo 
mass of the satellite is then set equal to the mass of the 
halo in which it lived at redshift Zform)- The position of a 
halo in this plane, relative to the redshift-dependent filter- 
ing mass (shown as a solid curve), determines the degree 
of PhS it experiences. We examine faint galaxies in two lu- 
minosity ranges, — f3 < Mb — 5 log ft < —12 (dots) and 
— 15 < Mb — 5 log ft < — 14 (small squares). 

Consider first a model with no PhS (left-hand panels). 
In this case, faint galaxies currently hosted in a 10 14 ft -1 Mq 
cluster halo (lower panel) have a much more extended range 
of formation redshifts than galaxies in a lO 12 ft _1 M Milky 
Way- type halo (upper panel). The fainter galaxies located 
in lower mass halos have a more extended range of forma- 
tion redshifts than their brighter counterparts, although this 
difference is rather small. The values of a for the two pop- 
ulations, given in each panel, indicate steep slopes in both 
environments. The filtering mass has no influence on these 
results, of course, but is shown to indicate which of ha- 
los might be expected to still form a galaxy when PhS is 
switched on. 

In the right-hand panels we plot results for a model with 
PhS. The effect of the filtering mass can now be clearly seen: 
halos below the solid line form galaxies much less efficiently 
than before (although there are certainly still some halos 
that managed to make galaxies below this line). The values 
of a given in the panels indicate that some flattening of the 
luminosity function has occurred, but this effect is small. 

Let us now examine the importance of these various 
processes using simplified models of PhS. Figure A2 shows 
luminosity functions in different mass halos constructed us- 
ing the simplified models described below, together with our 
standard model of PhS (heavy solid line) and a model with 
no PhS (thin solid line). 

Model A In the simplest scenario, suppose the filtering 
mass jumps abruptly from zero to Mf at redshift z r cion, 
and that any halo of smaller mass forming at lower redshift 
makes no galaxy at all, while galaxy formation in other ha- 
los proceeds unchanged. In this case, the luminosity func- 
tion should be unaffected at bright magnitudes (since these 
galaxies form in halos which are more massive than Mf), 
should plummet sharply at a magnitude corresponding to 
halo mass equal to Mf, and should rise again to fainter 
magnitudes (since these galaxies inhabit lower mass halos 
which typically form at redshifts greater than z ro ion)- On 
the basis of eqn. (A2), the slope of this rise should be large 
for low Mhost halos since these have the greatest variation in 
progenitor number with mass at z rc i on , and should approach 
some fixed value for large Mhost- This behaviour is borne out 
when we apply this simple prescription to our semi-analytic 
models, using values of Mf = 3 x 10 10 /i _1 Mq (compara- 



ble to the value of Mf at the redshifts where the major- 
ity of the faint galaxies are formed in our standard model) 
and z ro ion = 6.5 (the redshift of reionization in our stan- 
dard model; dot-long-dashed lines in Fig. A2) — suppression 
is greatest for small Mhost, tends to steepen the faint end 
of the luminosity function, and causes the most steepening 
for the lowest Mhost- This model produces a discontinuity in 
the luminosity function and is a poor approximation to the 
full treatment. 

Model B At the next level of sophistication, we take into 
account the variation of the filtering mass with redshift as 
predicted by our model, but retain the simple prescription 
wherein any halo forming with mass below the filtering mass 
forms no galaxy at all. Since Mf increases with time, galaxy 
formation in lower mass halos is suppressed earlier, an ef- 
fect which should tend to flatten luminosity function slopes. 
Since halos of a given mass form earlier in larger Mhost, 
we still expect the effects of filtering to be weaker for clus- 
ters than for the Local Group. At a given magnitude, this 
process would simply reduce the number of galaxies. How- 
ever, the detailed consequences for the faint end slope are 
now harder to anticipate since they depend on just how Mf 
varies with redshift. Applying this prescription to our model 
(long-dashed lines in Fig. A2) , we still find a discontinuity in 
the luminosity function (although this is much weaker than 
in model (i) since Mf varies smoothly with time) which is 
flattened relative to that without any PhS. Furthermore, the 
slopes are significantly flatter than in case (i), indicating that 
the time variation of Mf largely counteracts the steepening 
of the luminosity function produced by the variation in halo 
formation distributions. No environmental variation is intro- 
duced by these processes. The resulting luminosity functions 
are quite similar to those from our full PhS calculations, al- 
though the level of suppression is somewhat greater. Again, 
the effects are largest for the lowest Mhost • 

Model C Alternatively, we can consider a prescription 
which has no time variation in Mf, but in which, instead, 
the degree of suppression varies smoothly with halo mass (as 
prescribed by Gnedin 2000)T We find that this model also 
produces luminosity functions similar to those from our full 
PhS calculations (dot-short-dashed lines), and with about 
the right degree of suppression for suitable choices of Mf 

and Zrcion ■ 

Model D We next consider a model in which Mf varies 
with time and in which the degree of suppression varies 
smoothly with halo mass. This then differs from our full 
PhS calculations only in neglecting suppression due to pho- 
toheating by the ionizing background. In this model, pho- 
toionization redistributes galaxies which previously had a 
given magnitude over a range of fainter magnitudes. At a 
given magnitude, some galaxies are lost due to suppression, 
but also new ones appear as brighter galaxies are partially 

^ Somcrvillc (2002) considered a similar model, but kept the 
circular velocity corresponding to Mp constant, resulting in a 
filtering mass that increased with time. This prescription, also 
adopted by Bullock, Kravtsov & Weinberg (2000) provides a 
better match to the results of our standard model. Here, however, 
we are intersted only in exploring simplified models. 
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Figure Al. Positions of galaxies in the halo-mass/formation-rcdshift plane. Upper and lower panels show results for galaxies which 
reside in 10 12 h -1 Mq and 10 14 fo — *Mq halos at z = respectively. Left and right-hand panels show results for models ignoring PhS 
and including PhS respectively. The solid line shows the filtering mass as a function of redshift. Dots indicate galaxies in the range 
— 13 < Mb — 51og/i < — 12 while small squares indicate galaxies in the range —15 < Mb — 5 log h < —14. An estimate of a based on 
these two magnitude bins is shown in each panel. 



suppressed, thus becoming fainter. The net effect is there- 
fore difficult to judge since it depends on the time variation 
of the filtering mass and on the previous shape of the lumi- 
nosity function. We can, however, safely say that the effect 
is still expected to be larger for lower Mhost- Applying this 
prescription we reproduce the results of our full PhS calcu- 
lation rather well (dotted lines in Fig. A2), finding slightly 
less suppression overall than in case (ii), but still with no 
environmental variation. 



Model E A smoother variation of suppression with halo 
mass would result in less flattening of the faint end slope. 
For example, the short-dashed lines in Fig. A2 show results 
for a suppression of the form 



/ b M h 



1 + M F /M h ' 



(A3) 



with the same variation of filtering mass with redshift as in 
our standard model. 
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Figure A2. B-band galaxy luminosity functions in halos of different mass. Each panel shows the mean predicted model luminosity 
function in an ensemble of dark matter halos of mass given in each panel (ranging from halos containing small galaxies to rich clusters). 
All models include the effects of tidal limitation. The heavy solid lines show results from our complete model of PhS, while thin solid 
lines show a model with no PhS. Dot-long dashed lines have a filtering mass which jumps abruptly from zero to 3 X 10 10 h- 1 M Q at 
z = 6.5 and a sharp suppression at that mass at lower redshifts (Model A). Long-dashed lines have time varying filtering mass, but a 
sharp suppression at that mass (Model B). The dot-short dashed line shows a model with a fixed filtering mass but a smooth variation 
of suppression with halo mass (Model C). Dotted lines are a simplified model with time- varying filtering mass and a smooth variation 
of suppression with halo mass (Model D). Short-dashed lines are the same, but with an extra-smooth variation of suppression with halo 
mass (Model E). 
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In conclusion, a simple model in which galaxy formation 
at z < 2roion is entirely suppressed below some fixed mass 
scale grossly overestimates the effects of PhS, and produces 
very different luminosity functions from those predicted by 
our full calculation. Including either a smoothly varying de- 
gree of suppression or a time- varying filtering mass produces 
better agreement with our full calculations, indicating that 
these two ingredients are of comparable importance. The 
degree of suppression (as characterised by the reduction in 
amplitude of the LF) for a given time dependence in Mf de- 
pends on the functional form used to characterize the vari- 
ation of suppression with halo mass. The sharper the tran- 
sition from weak to strong suppression, the greater the net 
effect on the luminosity function. Smoother transitions tend 
to produce steeper faint end slopes, but even the sharpest 
transition possible does not produce slopes as flat as are 
observed and we never find any significant environmental 
variation. 
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